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EFFECTS  OF  SURFACE  REFLECTIONS 
OH  SHOCK  WAVE  IMPULSE 

by 

Verna  K.  Shuler 


ABSTRACT:  Under  free-field  conditions,  the  impulse  of  an  underwater  explosion  shock 
wave  is  constant  at  a  given  radial  distance  from  Idle  explosive  charge.  If  either 
the  charge  or  recording  gage  is  at  shallow  depth,  however,  the  rarefaction  wave 
reflected  back  into  the  vat.^r  from  the  interface  (the  "surface  cot. -off  wave") 
decreases  the  duration  of  the  positive  shock  vrve  so  that  the  radius  of  constant 
impulse  is  significantly  decreased  from  its  fTse-mter  value.  The  effect  of  surface 
cut-off  is  discussed  and  illustrated  in  this  report.  An  equation  for  shock  wave 
"partial  impulse,"  which  takes  into  account  the  effects  of  surface  cut-off,  is 
discus  pM  and  a  computer  program  that  generates  a  locus  of  points  of  constant  partial 
impulse  is  described. 
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The  impulse  of  an  underwater  explosion  shock  wave,  a  parameter  that  is  often  impor¬ 
tant  in  determining  the  damaging  effect  of  the  explosion,  may  be  reduced  to  a  small 
fraction  of  its  free -water  value  if  either  the  charge  or  the  target  is  at  very 
shallow  depth.  This  report  discusses  a  semi -empirical  equation  that  can  be  used  to 
approximate  the  near- surf ac»  shock  wave  impulse  for  various  test  configurations. 

An  iterative  machine  program  for  estimating  contours  of  constant  impulse  is  also 
presented.  This  program,  which  is  flexible  enough  to  be  used  for  a  number  of  appli¬ 
cations,  is  a  significant  improvement  over  the  tedious  and  time-consuming  manual 
methods  previously  used  for  estimating  such  contours. 
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EFFECTS  OF  SURFACE  REFLECTIONS  ON  SHOCK  HAVE  IMFUXfiE 


1.  INTRODUCTION 

From  various  underwater  explosion  tests  conducted  in  the  past,  empirical  rela¬ 
tionships  representing  the  shock  wave  impulse  for  particular  situations  have  been 
developed.  Arons  (references  a  and  b  ^describes  the  impulse  at  various  ranges  from 
a  shallow  charge  in  free  water  as  a  function  of  charge  weight  and  rcnge.  His  empiri¬ 
cal  power  functions  were  derived  from  pres sure -time  records  integrated  over  relatively 
long  times,  to  a  point  where  the  pressure  had  decayed  to  low  amplitude.  Slifko 
(reference  c)  has  developed  an  empirical  relationship  for  the  impulse  of  the  entire 
first  positive  phase  of  pressure  waves  measured  at  long  ranges  from  deep  charges  in 
free  water. 

A  situation  frequently  of  interest  in  practical  applications,  for  which  no 
provision  is  made  in  either  of  the  above-mentioned  relationships,  is  one  in  which 
surface  reflections  arrive  at  the  gage  shortly  after  the  shock  wave  arrival.  When 
the  negative  pressure  in  the  surface-reflected  wave  interacts  with  the  positive 
direct  pressure  wave,  it  cancels,  or  "cuts-off,  ”  the  pressure  in  the  tail  of  the 
shock  wave,  thereby  decreasing  the  impulse  at  the  gage  location. 

Because  surface  reflections  significantly  reduce  the  shock  wave  impulse  when 
either  the  charge  or  the  gage  is  relatively  shallow,  a  semi-empirical  equation  for 
shock  wave  impulse  that  includes  the  effects  of  surface  reflections  has  been  devel¬ 
oped  and  is  discussed  in  this  report.  This  quantity  is  called  "partial"  impulse  to 
differentiate  it  from  the  free  water  relationship  given  by  Arons.  Like  that  of 
Arons,  it  is  a  function  of  charge  weight  and  range,  hut  the  time  to  vixich  the  pulse 
is  integrated  is  included  as  a  variable  in  the  equation. 

In  addition  to  the  discussion  of  the  partial  impulse  equation,  the  report 
includes  a  computer  program  which  computes  the  locus  of  points  at  which  the  under¬ 
water  shock  wave  impulse  is  constant.  Instructions  for  using  the  program  are  also 
included. 


2.  SEMI-EMPIRICAL  EXPRESSION  FOR  SHOCK  WAVE  IMPULSE 
2.1  Assumed  Wave  Form 

The  pressure  in  the  shock  wave  from  an  underwater  explosion  rises  almost 
instantaneously  to  a  maximum  and  then  decays  to  ambient  pressure  at  a  rate  that 
depends  on  the  weight  of  the  explosive  charge  and  the  distance  the  wave  has  travelled 
outward.  The  initial  rate  of  decay  can  be  represented  quite  accurately  by  an  e  po- 
nential;  at  later  times,  however ,  the  pressure  decays  more  slowly  so  that  the  approxi¬ 
mation  of  the  wave  form  by  the  initial  exponential  underestimates  the  pressure  in 
the  tail  of  the  shock  wave.  Since  this  report  is  concerned  primarily  with  situations 
in  which  the  surface  cut-off  occurs  soon  after  the  shock  wave  arrival,  the  error 
introduced  by  using  the  exponential  approximation  is  negligible  for  cany  conditions 
of  practical  interest.  The  pressure  wave  for  which  impulse  is  computed  in  this 
report  is  assumed  to  be  an  exponential  of  the  form: 


References  are  given  on  page  iii. 


(1) 


P(t)  =  exp  (-  t/0) 

m 


itiere  the  peek  pressure,  Pm!  end  time  decay  constant,  9,  conform  to  the  fH?  simili¬ 
tude  values  given  fcy  Arons  (reference  a); 


P  (psi)  *  2.1 6  x  10k  (W1^3^)1*13 


(2) 


0  (msec)  =  6  x  10"2  W3  (W3/!?)"0*22  (3) 

vhere:  W  *  charge  weight  (lb) 

R  =  slant  range  (ft). 


2.2  Partial  Impulse 

The  shock  wave  impulse,  I,  is: 


1  *  y* 1  p  dt  (u) 

o 


and  from  Equations  1  and  4,  the  partial  Impulse,  I  ,  as  a  function  of  the  integration 
time,  is;  * 

\  *  ?m  3  l1  '  exp  ^l/e';] 


and  from  Equations  (2)  and  (3}: 


I  (psi-asec)  *  1300  W°*&  R-0*91  fl  -  exp  (Q) 
P  L  J 


vhere:  Q  =  -  16.67  tjW-0*26  R~0,22. 


(6) 


Integration  time,  t^,  for  the  partial  impulse  may  he  either  scene  specified  constant 
time,  t0,  or  surface  cut-off  time,  tg.  Since  the  surface  reflected  wave  of  negative 
amplitude  cancels  the  positive  pressure  in  the  shock  wave  as  shown  in  Figure  1, 
maximum  possible  integration  time  is  surface  cut-off  time. 

2.3  Surface  Cut-Off  Time 

Assuming  regular  acoustic  reflection  and  constant  wave  propagation  velocity, 
surface  c”toif  time,  tp,  is  determined  by  the  particular  charge -gage  geometry. 

From  the  geometric  relationships  shown  Figure  2,  the  following  equation  for  sur¬ 
face  cut-off  time  can  he  derived  if  the  velocity  of  sound  in  the  water  is  assumed  to 
he  5  ft /msec: 
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tg  (msec)  =  0.2  [(R2  +  4yD)1/2  -  e]  (7) 

•where:  R  =  slant  range  from  charge  to  gage,  ft 
Y  =  gage  depth,  ft 
D  =  charge  depth,  ft. 


An  equation  has  "been  derived  for  the  locus  of  gage  positions  at  vhich  surface  cut¬ 
off  time  is  constant.  Such  a  locus  is  the  lower  branch  of  a  hyperbola  with  a  verti¬ 
cal  transverse  axis  through  the  charge.  The  general  equation  for  the  hyperbola  is: 


(8) 


where:  a 
b 
X 


l/2  £R  p  ?  ,/? 

1/2  (to2  -  ar2)1'2 

horizontal  range 

difference  in  path  length  of  the 

the  direct  wave  (also  AR  =  t  V, 

water). 


reflected  wave  and  the  path  length  of 
where  V  s  velocity  of  sound  in  the 


Figure  3  shows  several  contours  of  constant  surface  cut-off  time  for  a  charge 
at  865-ft  depth.  At  the  water  surface  tc  -  0;  on  each  narrowing  hyperbolic  branch 
t  increases  and  finally  reaches  a  maximum  value  of  0.346  sec  when  the  gage  is 
directly  below  the  charge.  Since  tfl  is  determined  from  purely  geometric  considera¬ 
tions,  the  charge  and  gage  locations  may  be  used  interchangeably.  In  other  words, 
if  the  gage,  rather  than  the  charge,  is  held  at  a  constant  depth  of  865  feet.  Figure 
3  shows  how  surface  cut-off  time  vaides  as  charge  location  changes. 

2.4  Effect  of  Surface  Cut-Off  Time  on  Partial  Impulse 

Figures  4  and  5  illustrate  the  effect  of  surface  cut-off  on  shock  wave  impulse 
at  a  given  horizontal  distance  from  the  charge.  In  these  figures,  10- lb  charges  are 
located  at  2- ft  and  20- ft  depths,  respectively,  and  horizontal  range  is  held  at 
constant  20  feet.  Maximum  pressure  (Equation  2)  ana  partial  impulse  (Equation  6), 
each  normalized  to  its  value  at  1-ft  depth,  is  plotted  versus  gage  depth;  integration 
time  is  surface  cut-off  time. 

In  Figure  4,  as  the  gage  is  lowered  from  the  surface  down  to  about  6  feet,  peak 
pressure  changes  very  little  because  slant  range  varies  only  slightly.  Partial 
impulse,  however,  increases  significantly  with  increasing  gage  depth  until  it  reaches 
a  maximum  of  about  4.5  times  its  1-ft  value  when  the  gage  is  about  13  feet  deep. 
Increasing  surface  cut-off  time  is  responsible  for  the  increase  in  impulse  until  the 
gage  reaches  about  13  feet;  at  greater  depths,  slant  range  increases  enough  to  off¬ 
set  the  surface  cut-off  effects  and  impulse  decreases  sharply.  The  pressure  also 
drops  off  more  rapidly  for  the  greater  gage  depths  because  of  the  increasing  slant 
range. 
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With  the  deeper  chsrge  (20  feet)  of  Figure  5,  there  is  less  variation  in  shock 
wave  impulse  as  the  gage  is  lowered.  The  impulse  here  varies  in  the  same  fashion  as 
peak  pressure;  both  parameters  are  at  a  maximum  vhen  the  charge  and  gage  are  at  the 
same  20-ft  depth.  Impulse  is  not  greatly  affected  by  surface  cut-off  when  the 
charge  is  at  20  feet  because  surface  reflections  arrive  relatively  late. 


3.  CONTOURS  OF  CONSTANT  3MHJLSF. 

3.1  Modification  of  Free -Fie Id  Contour 

Under  free-field  conditions,  the  contours  of  constant  impulse  lie  along  spheri¬ 
cal  surfaces  centered  at  the  explosion.  For  a  given  weight  of  charge,  the  radius  s.t 
which  the  free-field  Impulse  has  some  particular  constant  value  can  be  found  from 
Arons'  empirical  equation  for  shock  waves  integrated  to  6. 7  6  (reference  b); 

I  (psi-msec)  =  1780  W1/3  (W^AO0-*  (9) 


Near  the  water  surface,  however,  the  radius  must  be  found  from  a  partial  impulse 
function  such  as  Equation  6.  At  some  point.  Equation  6  becomes  a  poor  approximation, 
as  discussed  below,  and  empirical  corrections  are  needed  to  relate  the  two  impulse 
functions.  Unfortunately,  only  fragmentary  data  obtained  at  relatively  short  ranges 
are  available  so  that  it  is  not  possible  to  empirically  define  contours  for  all 
conditions  of  interest.  A  sample  contour  -that  could  be  estimated  from  available 
data  is  shown  in  Figure  6. 

In  Figure  6,  the  heavy  line  is  the  estimated  contour  along  which  I  =  100 
psi-msec  for  a  10-Ib  charge  detonated  at  10  feet.  For  shallow  depths,  this  contour 
coincides  with  the  partial  impulse  contour  computed  from  Equation  6  of  this  report; 
at  greater  depths  the  contour  coincides  with  the  free-field  contour  of  Arons-  These 
two  limiting  curves  are  shown  dashed  beyond  the  point  where  they  coincide  with  the 
desired  contour.  Interpolation  between  the  two  limiting  contours  was  based  on  the 
fragmentary  data  mentioned  above. 

3-2  Comparison  of  Contours 

It  can  be  seen  from  Figure  6  that  for  these  particular  conditions,  the  partial 
impulse  equation  discussed  in  previous  sections  of  this  report  provides  a  poorer  fit 
as  the  gage  location  deepens,  and  gives  an  estimated  slant  range  that  is  about  15 $ 
too  small  when  the  gage  depth  is  UO  feet.  For  a  different  set  of  conditions, 
the  partial  impulse  equation  given  here  might  be  an  excellent  aT  roximation  down  to 
much  greater  depths,  but  there  is  no  simple,  straightforward  way  to  extrapolate 
Figure  6  to  other  conditions.  The  adequacy  of  the  exponential  approximation  must 
he  reckoned  in  terms  of  the  number  of  time  constants  (8  of  Equation  3)  over  tdiich 
the  wave  is  integrated,  and  there  is  no  simple  generalised  expression  for  this 
relative  quantity,  t  =  tj/0.  The  time  constant,  8,  depends  upon  only  chsrge  weight 
and  range,  while  t  depends  upon  charge  weight,  total  geometry,  level  of  constant 
impulse  of  interest  and  upon  the  maximum  integration  time  if  such  a  limit  is  imposed. 
We  do  not  have  adequate  data  to  define  empirical  corrections  to  Equation  6  for  all 
cases  of  interest;  however,  certain  generalizations  can  be  made.  If  t,  is  no  greater 
than  about  3  0,  Equation  6  probably  gives  an  adequate  approximation  ofl^.  In  terms 
of  the  variables  noted  on  Figure  6,  this  means  that  the  exponential  approximation 
will  fit  the  true  contour  down  to  deeper  gage  depths  if  (a)  charge  weight  is  increased 
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or  (b)  the  impulse  level  of  the  desired  contour  is  decreased  or  (c)  the  charge  is 
fired  at  a  shallower  depth,  When  sufficient  data  are  available,  generalized  methods 
of  generating  the  contours  of  constant  impulse  for  a  range  of  variables  will  be 
developed. 


4.  COMPUTER  PROGRAM  FOR  ESTIMATING  LOCI  OF  COBST WfT  PARTIAL  IMPULSE 

4.1  General 

Since  the  range  corresponding  to  a  constant  partial  impulse  cannot  be  found 
explicitly  from  Equation  6,  estimating  loci  of  constant  partial  impulse  by  manual 
computation  is  tedious  and  time-consuming.  With  the  iterative  computer  program 
outlined  below,  however,  the  range  at  *Aich  I  is  within  ±0.1  psi-msec  of  some 
constant  value,  IQ,  can  be  found  rapidly.  p 

4.2  Outline  of  Program 

The  program  has  been  written  in  the  Basic  Language  on  the  remote  computing 
facility  at  SOL.  It  finds  the  maximum  range  to  which  a  charge  of  weight  W,  at  depth 
D,  can  deliver  some  given  Impulse  IQ,  within  some  specified  tins  tQ.  Although 
tQ  is  introduced  as  a  constant  value  in  the  program  listing  shown  here,  it  can  be 
introduced  in  functional  form  with  only  minor  program  changes.  Note  that  the  maxi¬ 
mum  possible  duration  of  the  shock  wave  is  the  surface  cut-off  time,  tg,  of  Equation 
7;  consequently,  the  program  compares  the  input  parameter  tQ  with  t  for  the  geometry 
in  question  and  uses  t  rather  than  tQ  in  computing  the  partial  impulse  'Aenever 
tQ  >  tg.  If  the  desired  integration  time  is  always  t  ,  then  a  large  value  for  tQ 
(such  as  10  sec  )  should  be  inserted  in  the  appropriate  data  statement . 

A  detailed  explanation  of  the  program  logic  is  not  included  in  the  report 
because  sufficient  comments  are  scattered  throughout  the  program  so  that  one  who 
understands  the  programming  language  should  be  able  to  follow  the  logic .  The  input 
parameters  from  which  the  computer  completes  all  computations  are  as  follows: 

(1)  charge  weight,  W,  lb 

(2)  charge  depth,  D,  ft 

(3)  gage  depth,  Y,  ft 

(4)  given  Impulse,  IQ,  psi-msec 

(?)  given  integration  time,  tQ,  msec. 

Hie  iteration  involves  these  four  basic  steps: 

(1)  A  slant  range,  R,  is  automatically  selected  and  surface  cut-off  time,  tg, 
for  the  particular  geometric  situation  is  calculated. 

(2)  Partial  impulse,  I.  is  computed  using  this  same  slant  range,  R,  gage 
depth,  Y,  and  charge  weight,  W.  Integration  time,  t,  is  the  smaller  of 
t  or  tg. 

(3)  Partial  impulse  from  step  (2)  is  compared  with  the  given  impulse,  IQ; 
if  Ip  <  (l0  -  O.l),  then  range  is  automatically  decreased  or  If  Ip  > 

(l0  *  0.1),  range  is  automatically  increased.  Hie  new  range  is  W  and 


If,  for  seme  reason,  there  is  no  convergence  after  50  iterations,  a  message  indicating 
where  adjustments  should  be  made  will  be  printed.  If  the  conditions  are  such  that 
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there  is  no  solution,  this  also  will  he  indicated.  A  list  of  the  program  and  samples 
of  output  are  included  in  Appendices  A  and  B  to  this  report. 

•i  3  Data  Arrangement 

All  data  for  the  Basic  program  is  entered  in  data  statements  between  statements 
5000  and  7ACO.  the  statement  mothers  used  here  as  examples  are  chosen  for  illustra¬ 
tion  purposes;  the  user  say  number  his  data  statements  as  he  chooses: 

5001  Data  for  this  statement  should  be : 

(a)  number  of  charge  weights 

(b)  number  of  gage  depths 

(c)  charge  depth  (ft). 

5002  This  data  statement  contains  the  individual  charge  weights  (lb)  given — 
there  may  be  20  of  these  values;  the  number  of  charge  weights  is  indicated 
in  5001  (a). 

5003  Each  individual  gage  depth  (ft)  is  listed  in  this  statement — again,  there 
may  be  20  gage  depths  given.  The  number  of  gage  depths  given  is  specified 
in  5001  (b). 

5004  In  this  statement,  the  desired  impulse  (Iq,  psi-nsec)  and  integration 
time  (t0,  msec)  are  given,  respectively. 

The  data  arrangement  is  illustrated  in  the  following  example:  For  10- lb  and  100- lb 
charges  fired  at  UO-ft  depth,  determine  the  ranges  at  which  IQ  *  150  psi-msec  viien 
t  *  2.5  ssec.  Use  gage  depths  of  10,  20,  and  50  feet.  Appropriate  data  statements 
follow: 

5001  Data  2,  3,  40 

5002  Efeta  10,  100 

5003  Data  10,  20,  50 

5004  Data  150,  2-5 

Note  that  a  new  run  is  required  if  charge  depth,  given  Impulse,  or  given  integration 
time  is  changed. 
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t - 1 — i — rum - 1 - 1 — i — i  i  ftt; 

10  IB  CHARGE  AT  2  FT  DEPTH 
CONSTANT  HORIZONTAL  RANGE  =  20  FT 


PARAMETER  AT  PEP' 
PARAMETER  AT  1  FT 


1  lit* 


uBim  or  oBMsvm  mom* 


10  DIM  R( 50) >  1(50);  T(50);  Y(20);  W(20> 

50  READ  C;  G»  D 

55  REM  C=NUMBER  0F  CHG  WTS  G=N UMBER  OF  GAGE  DPHS  D=CKG  DPH 
60  FOR  L  =  1  T0  C 
70  READ  W(L) 

75  REM  W  =  INDIVIDUAL  CHARGE  WEIGHT;, 

80  NEXT  L 

90  FOR  M  =  1  T0  G 

100  READ  Y(M) 

105  REM  Y  =  INDIVIDUAL  GAGE  DEPTHS 
110  NEXT  M 
1  12  PRINT 
113  PRINT 

115  PRINT  "CHARGE  DEPTH  =  "  D"  FEET" 

118  READ  10; TO 

119  REM  10  =  GIVEN  IMPULSE  TO  =  GIVEN  INTEGRATION  TIME 

120  FOR  K  =  1  TO  C 
122  PRINT 

124  PRINT 

130  PRINT  "CHARGE  WEIGHT  (LB)  ="  W(K) 

140  FOR  M  =  1  TO  G 

141  LET  XI  =  1000 

143  LET  F  =  0 

144  LET  H  =  0 

145  PRINT 
147  PRINT 

150  PRINT"GAGE  DEPTH  (FT)  =  "Y(M) 

190  LET  T(M)  =  TO 
195  LET  A1  =  10+.! 

197  LET  B1  =  IO-.l 

200  PRINT  "GIVEN  IMPULSE  (PSI-MSEC)  =  "  10 
210  PRINT  "GIVEN  INTEGRATION  TIME  (MSEC)  =  "T(M) 

220  LET  J  =  1 

230  LET  R(J)  =  (4*Y(M)*D  -  (5*T0)t2>  /  (10  *  TO) 

233  IF  R(J)  <  0  THEN  400 

240  REM  RANGE  WHERE  SURFACE  CUT-OFF  TIME  =  GIVEN  INTEG.  TIME  IS  R(J) 

250  GOSUB  1000 

260  IF  I(J)  >=I0  THEN  500 

262  LET  XI  =  R(J) 

265  REM  MAX  RANGE  SO  FAR  IS  XI. 

270  REM  COMPUTED  IMPULSE  IS  LESS  THAN  GIVEN  IMPULSE. 

278  LET  F  =  1 
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280  F0R  J  =  2  T0  50 

290  LET  Q  =  J-l 

300  LET  RCJ)  =  • 5*RCQ) 

31 0  REM  TRYING  A  RANGE  1/2  0F  PREVIOUS  RANGE. 

320  G0SUB  1000 

330  IF  ICJ)  >=  10  THEN  700 

350  REM  FINDING  A  RANGE  SMALL  ENOUGH  S  'HAT  GIVEN  IMP  CAN  BE  DELIVERED. 
360  NEXT  J 

370  PRINT'*AFTER  50  ITERATIONS  0N  DECREASED  RANGE*  GIVEN  IMP  NOT  REACHED’* 

390  G0  T0  5000 

400  LET  K  =  1 

402  G0SUB  3000 

404  IF  Y(M)  =  D  THEN  426 

406  REM  CASE  WHERE  CHARGE  AND  GAGE  LIE  0N  A  VERTICAL  LINE. 

408  LET  RCJ)  =  ABSCYCM)  -  D) 

410  G0SUB  1000 

412  IF  ICJ)  >=  10  THEN  420 

415  PRINT”SINCE  IMPULSE  CALCULATED  AT  THE  LARGER  0F  GIVEN  INTEGRATION” 

416  PRINT”TIME  0R  SURFACE  CUT-OFF  TIME  IS  LESS  THAN  THE  GIVEN  DESIRED” 

417  PR  I  NT”  IMPULSE*  THERE  IS  N0  SOLUTION.  ” 

418  GO  TO  995 

420  REM  SINCE  COMPUTED  IMP>GIVEN  IMP.*  RANGE  MAY  BE  INCREASED. 

422  LET  N1  =  RCJ) 

424  REM  MIN  RANGE  SO  FAR  IS  Nl. 

425  GO  TO  427 

426  LET  Nl  =  0 

427  LET  XI  =  10 

428  IF  XI  «=N1  THEN  432 
430  GO  TO  434 

432  LET  XI  =  Nl  +  10 

434  LET  X9  =  1 

436  LET  X8  =  XI 

438  FOR  J  =  2  TO  50 

440  LET  RCJ)  =  X9  *  X8 

442  GOSUB  2000 

444  GOSUB  4000 

460  GOSUB  1000 

470  IF  ICJ)  >=  10  THEN  700 

472  LET  X8  =  X8  /  2 

474  NEXT  J 

480  PRINT”N0  CONVERGENCE  AFTER  50  ITERATIONS  ON  438»" 

485  GO  TO  5000 

500  REM  SINCE  COMPUTED  IMP>GIVEN  IMP  WHERE  SURFACE  CUT-OFF  TIME  EQUALS 
505  RF.M  GIVEN  INTEGRATION  TIME*  ATTEMPT  IS  MADE  TO  FIND  MAX  RANGE  AT 
510  REM  WHICH  GIVEN  IMP  CAN  BE  DELIVERED  WITHIN  SURF  CUT-OFF  TIME. 

534  LET  Nl  =  RCJ) 

535  REM  MIN  RANGE  IS  Nl . 

540  LET  U  =  10 

545  LET  Z3  =  RCJ) 

550  FOR  J  =  2  TO  50 

570  REM  TRYING  A  RANGE  OF  U  *  Z3. 

580  LET  RCJ)  =  U*Z3 
590  GOSUB  2000 
600  GOSUB  1000 
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olO  IF  ICJ)>=I0  THEN  700 
617  LET  XI  =  RCJ) 

619  REM  MAX  RANGE  IS  XI. 

620  REM  COMPUTED  IMPULSE  LESS  THAN  GIVEN  IMPULSE. 

630  LET  U  =  U/2 

640  NEXT  J 

650  PRINT  "NO  CONVERGENCE  AFTER  50  ITERATIONS  ON  PARAG.  550** 

660  GO  TO  5000 
700  LET  N1  =  RCJ) 

710  REM  MIN  RANGE  SO  FAR  IS  Nl. 

712  IF  XI  «=  Nl  THEN  716 

71 4  GO  TO  720 

716  LET  XI  =  Nl  *100 

718  REM  LETTING  MAX  RANGE  BE  XI. 

720  FOR  J  =  3  TO  50 

730  LET  RCJ)  =  XI  -  CXl-ND/2 

740  REM  TRYING  A  RANGE  1/2  OF  WAY  BETWEEN  Nl  AND  XI. 

750  IF  F  =  1  THEN  780 
760  GOSUB  2000 
770  GO  TO  790 
780  LET  TCM)  =  TO 
790  IF  H< 1  THEN  799 
792  GOSUB  4000 

799  GOSUB  1000 

800  IF  ICJ)  >  B 1  THEN  825 
810  GO  TO  835 

825  IF  ICJ)  «  A1  THEN  895 
835  IF  ICJ)  >  10  THEN  860 

837  REM  COMPUTED  IMP  NOT  WITHIN  GIVEN  IMP  +0R-  0.1  IMP  TOO  SMALL. 

845  LET  XI  =  RCJ) 

850  REM  NEW  MAX  RANGE  IS  XI. 

855  GO  TO  880 

860  REM  COMPUTED  IMP  NOT  WITHIN  GIVEN  IMP  +  0R-  0.1  IMP  TOO  LARGE. 

865  LET  Nl  =  RCJ) 

870  REM  NEW  MIN  RANGE  IS  Nl. 

880  NEXT  J 

885  PRINT  "NO  CONVERGENCE  AFTER  50  ITERATIONS  ON  720." 

890  GO  TO  5000 

895  LET.  K4  ~  D- Y  CM) 

896  IF  RCJ)  <ABSCK4)  THEN  985 

980  IF  H  =  1  THEN  983 

981  IF  F  =  1  THEN  990 

982  GO  TO  987 

983  IF  TCM)  <  TO  THEN  987 

984  GO  TO  990 

985  PRINT"N0  SOLUTION - GAGE  AND  CHARGE  ARE  TOO  FAR  APART." 

986  GO  TO  995 

987  PRINT-RANGE  OF"RCJ)"FEET  HAS  BEEN  FOUND  AT  WHICH  IMPULSE  0F"ICJ) 

988  PR I NT" CAN  BE  DELIVERED  IN"TCM)"MSEC  WHICH  IS  SURFACE  CUT-OFF  TIME. 

989  GO  TO  995 

990  PRINT"RANGE  OF"RCJ)"FEET  HAS  BEEN  FOUND  AT  WHICH  IMPULSE  0F"ICJ) 

99 1  PRINT"CAN  BE  DELIVERED  WI THIN"TCM)"MSEC  WHICH  IS  GIVEN  INTEGRATION 

992  PRINT  "TIME.  SURFACE  REFLECTIONS  DO  NOT  ARRIVE  THIS  EARLY." 

99  5  NEXT  M 
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996  NEXT  K 
998  G0  T0  7000 

1000  REM  THIS  SUBROUTINE  COMPUTES  IMPULSE. 

1030  LET  0=-16.67*TCM)/WCK)t.333*  CWCK)*.333/  RtJ))*.22 

1040  LET  ICJ)  =  1 300* WOO  t . 333*C  WCK)  t . 333/  RC J) ) t .9 1 *C 1-EXPCQ) > 

1090  ’RETURN 

2000  REM  THIS  SUBROUTINE  COMPUTES  SURFACE  CUT-OFF  TIME  FOR  EACH  RANGE. 
2030  LET  TCM)=CCRCJ)t2+4*YCM)*D)t.5-RCJ))  /  5 
2060  RETURN 

3000  REM  THIS  SUBROUTINE  COMPUTES  MAX  TS  C SURFACE  CUT-OFF  TIME). 

3001  REM  CTO  IS  GIVEN  INTEGRATION  TIME.) 

3002  REM  NO  SOLUTION  FOR  TS  =  TO  MUST  FIND  MAX  TS  AND  PROCEED. 

300  5  IF  D  <  YCM)  THEN  3020 

3010  LET  TCM)  =  2*YCM)  /  5.0 

3011  IF  D  =  YCM)  THEN  3017 
3015  GO  TO  3035 

3017  RETURN 

3020  LET  TCM)  =  2*D  /  5.0 
3035  IF  TCM)  <  TO  THEN  3045 
3040  LET  TCM)  =  TO 

3042  REM  TO<MAX  TS  USE  TO  IN  COMPUTING  IMPULSE* 

3043  RETURN 

3045  REM  MAX  TS<TO  USE  TS  IN  COMPUTING  IMPULSE. 

3050  RETURN 

4000  REM  THIS  SUBROUTINE  FINDS  SMALLER  OF  TO  AND  TS  TO  BE  USED  FOR  IMP” 
4010  IF  TO  >  TCM)  THEN  4050 
402  0  LET  TCM)  -  TO 

4030  REM  SINCE  T0<TS*  TO  WILL  BE  USED  IN  CALCULATING  IMPULSE. 

4040  RETURN 

4050  REM  SINCE  TS  <  TO*  TS  WILL  BE  USED  IN  CALCULATING  IMPULSE. 

4060  RETURN 

5000  PRINT  "ITERATIONS  NEEET  REVISING" 

5001  DATA  2  3*40 

5002  DATA  10*100 

5003  DATA  10*20*50 

5004  DATA  150*2.5 
7000  STOP 

8000  END 
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CHARGE  WEIGHT  CLB3  =  10 


GAGE  DEPTH  CFT3  =  10 

GIVEN  IMPULSE  EPSI-MSEC3  =  150 

GIVEN  INTEGRATION  TIME  CMSEC3  =  2.5 

RANGE  OF  53.6331  FEET  HAS  BEEN  FOUND  AT  WHICH  IMPULSE  0F  150.023 
CAN  BE  DELIVERED  WITHIN  2.5  MSEC  WHICH  IS  GIVEN  INTEGRATION 

TIME.  SURFACE  REFLECTIONS  DO  NOT  ARRIVE  THIS  EARLY. 


GAGE  DEPTH  CFT3  =  20 

GIVEN  IMPULSE  CPSI-MSEC3  =  150 

GIVEN  INTEGRATION  TIME  CMSEC3  =  2.5 

RANGE  OF  53.6223  FEET  HAS  BEEN  FOUND  AT  WHICH  IMPULSE  OF  150.051 
CAN  BE  DELIVERED  WITHIN  2.5  MSEC  WHICH  IS  GIVEN  INTEGRATION 

TIME.  SURFACE  REFLECTIONS  DO  NOT  ARRIVE  THIS  EARLY. 


GAGE  DEPTH  CFT3  =  50 

GIVEN  IMFJLSE  C PSI-MSEC3  =  150 

GIVEN  INTEGRATION  TIME  CMSEC3  =  2.5 

RANGE  OF  53.629  FEET  HAS  BEEN  FOUND  AT  WHICH  IMPULSE  OF  150*034 
CAN  BE  DELIVERED  WITHIN  2.5  MSEC  WHICH  IS  GIVEN  INTEGRATION 

TIME.  SURFACE  REFLECTIONS  DO  NOT  ARRIVE  THIS  EARLY. 


CHARGE  WE*  JHT  CLB3  =  100 


GAGE  DEPTH  CFT3  =  10 

GIVEN  IMPULSE  E  PSI-MSEC3  =  150 

GIVEN  INTEGRATION  TIME  CMSEC3  =  2.5 

RANGE  OF  190.62  FEET  HAS  BEEN  FOUND  AT  WHICH  IMPULSE  OF  149.909 
CAN  BE  DELIVERED  IN  .830324  MSEC  WHICH  IS  SURFACE  CUT-OFF  TIME. 
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GAGE  DEPTH  CFT3  =  20 

GIVEN  IMPULSE  CPSI-MSEC3  =  ISO 

GIVEN  INTEGRATION  TIME  CMSEC3  =  2.5 

RANGE  0F  231.402  FEET  HAS  BEEN  FOUND  AT  WHICH  IMPULSE  0F  149.994 
CAN  BE  DELIVERED  IN  1.36281  MSEC  WHICH  IS  SURFACE  CUT-OFF  TIME. 


GAGE  DEPTH  CFT3  =  50 

GIVEN  IMPULSE  CPSI-MSEC3  =  150 

GIVEN  INTEGRATION  TIME  CMSEC3  =  2.5 

RANGE  0F  261.05  FEET  HAS  BEEN  F0UND  AT  WHICH  IMPULSE  0F  149.955 
CAN  BE  DELIVERED  WITHIN  2.5  MSEC  WHICH  IS  GIVEN  INTEGRATION 
TIME.  SURFACE  REFLECTI0NS  D0  N0T  ARRIVE  THIS  EARLY. 
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